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Abstract 

We present a derivation of the threshold resummation formula for the Drell-Yan rapidity 
distribution. Our argument is valid for all values of rapidity and to all orders in perturbative 
QCD and can be applied to all Drell-Yan processes in a universal way, i.e. both for the production 
of a virtual photon 7* and the production of a vector boson W^, . We show that for the fixed- 
target experiment E866/NuSea used in current parton fits, the NLL resummation corrections 
are comparable to NLO fixed-order corrections and are crucial to obtain agreement with the 
data. 
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In perturbative QCD, it is well known that, when one approaches to the boundary of the phase 
space, the cross section receives logarithmically-enhanced contributions at all orders. These 
large terms have been resummed a long time ago for the classes of inclusive hadronic processes 
of the type of deep- inelastic and Drell-Yan IH [2] to next-to-leading- logarithmic order (NLL). 
More recently, the next-to-next-to-leading-logarithmic (NNLL) accuracy has been reached [3]. 
Threshold resummation of inclusive processes can affect significantly cross sections and the 
extraction of parton densities [H |5] . For the case of small transverse momentum distributions in 
Drell-Yan processes, it has been shown that resummation is necessary to reproduce the correct 
behavior of the cross section [6]. 

The differential rapidity Drell-Yan cross section is used for the extraction of the ratio d/u 
of parton densities. The accurate knowledge of these functions is needed to study Higgs boson 
production and the asymmetry . The resummation of Drell-Yan rapidity distributions was 
first considered in 1992 [7]. At that time, it was suggested a resummation formula for the case of 
zero rapidity. Very recently, thanks to the analysis of the full NLO calculation of the Drell-Yan 
rapidity distribution, it has been shown [8], that the result given in [7j is valid at NLL for all 
rapidities. 

In this Letter, we will give a simple proof of an all-order resummation formula valid for all 
values of rapidity. To do this, we will use the technique of the double Fourier-Mellin moments 
developed in [9|. In particular, we will show that the resummation can be reduced to that of the 
rapidity-integrated process, which is given in terms of a dimensionless universal function for both 
DY and and production, and has been largely studied [U [2] even to all logarithmic orders 
jlOj . Finally, we implement numerically the resummation formula and give predictions of the 
full rapidity-dependent NLL Drell-Yan cross section for the case of the fixed-target E866/NuSea 
experiment. We find that resummation at the NLL level is necessary and that its agreement 
with the experimental data is better than the NNLO calculation 

We consider the general Drell-Yan process in which the collisions of two hadrons {Hi and 
H2) produce a virtual photon 7* (or an on-shell vector boson V) and any collection of hadrons 



In particular, we are interested in the differential cross section ^^^(a;, Q^, Y)^ where Q"^ is the 
invariant mass of the photon or of the vector boson, x is defined as usual as the fraction of 
invariant mass that the hadrons transfer to the photon (or to the vector boson) and Y is the 
rapidity of 7* [V) in the hadronic center-of-mass: 



where E andp^ are the energy and the momentum along the collisional axis of 7*(V") respectively. 
At the partonic level, a parton 1(2) in the hadron Hi {H2) carries a longitudinal momentum 
Pi = XiPi (p2 = 2:2^2) • Thus, the rapidity in the partonic center-of-mass (y) is obtained 
performing a boost of Y between the two frames: 



(X): 



HiiPi) + H2{P2) ^ 1*{V){Q) + X{K). 



(1) 





(3) 
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In order to understand the kinematic configurations in terms of rapidity, it is convenient to 
define a new variable u, 



Q ' Pi -2v -^1 -2Y r A \ 

u = — = e ^ = — e . (4) 

Q ■P2 X2 



which can assume all the values in the closed interval, 



1 

z<u<-, (5) 
z 



with 



2pi-J92 (Pl+P2)^ 2:1X2' 



(6) 



The upper and lower bounds in eq.([5]) are reached when the extra radiation is emitted collinear 
to the incoming parton 1 and 2 respectively. Eqs. (!3|^ allow us to rewrite the relation in eq.([5]) 
as a relation for the upper and lower bounds of the partonic center-of-mass rapidity: 

1 , 11 

- < y < - In -. (7) 

Substituting eqs.(jll[6]) into the two conditions u > z and u < 1/z, we obtain the lower kinematic 
bound for Xi and X2'- 



Xi > \fxe^ = Xi, X2 > y/xi 



-Y 



xe = X2 [o) 

and the obvious requirement x'j'(2) ^ 1 implies that the hadronic rapidity has a lower and an 
upper bound: 

-lnx<r<-ln-. (9) 

2 - - 2 X ^ ' 

The variable z in eq. can be viewed as the fraction of invariant mass that the incoming partons 
transfer to 'j*{V). and, hence, the threshold limit is reached when z approaches to 1. 

According to the standard factorization of collinear singularities of perturbative QCD, the 
expression for the hadronic differential cross section in rapidity has the form. 



where y depends on Y, xi and X2 according to eq.Q. The sum runs over all possible partonic 
subprocesses, f/^-* ,Fj'^^ are respectively the parton densities of the hadron Hi and H2, /x is 
the factorization scale (chosen equal to renormalization scale for simplicity) and ddij/ {dQ'^dy) 
is the partonic cross section. In the threshold limit the gluon-quark channels are suppressed 
by powers of (1 — z) [10] and, so, in order to study resummation, we will consider only the 
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quark-anti-quark contributions of the sum in eq. fllOl) . These last terms are related to the same 
dimensionless coefficient function C(z, ^^(/i^), through the relations 

for the virtual photon vertex and 

^^^^ dqZy I ^i'^2,^,a.(/i ),yj = ^ S{Q - M^jC I 2;, — , ), ?/ 1 , (12) 

for the real vector boson vertex. Here Gi? is the Fermi constant, My is the mass of the produced 
vector boson. The coefficients Cgq' are given by: 

Cqq' = QlSgg for 7*, (13) 
c,q, = I VP forl^±, (14) 

Cqq' = mg'.r + i9ir]Sqq forZ", (15) 

where is the square charge of the quark q, Vqqi are the CKM mixing factors for the quark 
flavors g, q' and 

14 1 

gl = - - - sin^Ow, 9l = ^ for an up-type quark, (16) 

12 1 
= -2 + 3sin^^iy, fi'a = "2 for a down-type quark, (17) 

with 6w the Weinberg weak mixing angle. Thus, we are left with a dimensionless cross section 
of the form: 

Jxl Xi JxO X2 \ fl^ J 

where Fi and F2 are quark or anti-quark parton densities in the hadron Hi and H2 respectively. 
This shows the universality of resummation in Drell-Yan processes in the sense that only the 
quantity defined in eq.( fT8|) has to be resummed. 

We shall now show that the resummed expression of eq. (fT8|) is obtained by simply replacing 
the coefficient function C{z,Q'^/ fi^,as{^'^),y) with its integral over y, resummed to the desired 
logarithmic accuracy. To show this, we recall that resummation is usually performed in the 
space of the variable N, which is the Mellin conjugate of x, since Mellin transformation turns 
convolution products into ordinary products. In the case of the rapidity distribution, however, 
this is not sufficient. In fact, we see that the Mellin transform with respect to x, 

a{N, Y) = dxx^~^a{x, Q^ Y), (19) 
Jo 



does not diagonalize the double integral in eq. flTSl) . because the partonic center-of-mass rapidity 
y depends on xi and X2 through eq.(l3]). The ordinary product in Mellin space can be recovered 
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performing the Mellin transform with respect to x of the Fourier transform with respect to Y . 
Using eqs. fl9|7|) and the fact that the coefficient function must be symmetric in y, we find 



a(Ar,Q^M) = ^dxx^-' r^"^'^ dYe'^^^a{x,Q\Y) (20) 

"'0 "'log 

= Ff ^ {N + zM/2, /x2)Ff ^ (AT - zM/2, ^i^)C (^iV, ^, as{^?) , A^j , (21) 

where 

Fl'^iN ±tM/2,i^^) = dxx^-'^'^'^^Fl'^ix^fi^), (22) 

Jo 

cfiV,^,«,(/i2),Mj = 2 dzz""-' dycos{My)Ciz,^,as{fi'),yj. (23) 



The dependence on M, the Fourier conjugate of the rapidity y, originates from the parton 
densities, that depend on N ±iM/2, and from the factor of cos(My) in the integrand of eg. (1231) . 
This last dependence, however, is irrelevant in the large-A^ limit. Indeed, one can expand 
cos (My) in powers of y, 



cos{My) = 1 ^ + 0{MY)- (24) 

and observe that the first term of this expansion leads to a convergent integral (the rapidity- 
integrated cross section), while the following terms are suppressed by powers of [1 — z), since 
the upper integration bound is 

lni= = l(l-.) + 0((l-.)2). (25) 



Hence, up to terms suppressed by factors 1/A^, eq.( l23|) is equal to the Mellin transform of 
the rapidity-integrated Drell-Yan coefficient function that we call Ci{N,Q'^/ fi'^,as{fi'^)). This 
completes our proof. We get 

cr'^^^(iV,Q^M) = Ffi(iV + 2M/2,/i2)Ff2(iV-2M/2,/i2)Cr (^iV, ^, ^.(/x^)^ . (26) 

This is the main theoretical result of our Letter: it shows that, near threshold, the Mellin- 
Fourier transform of the coefficient function does not depend on the Fourier moments and that 
this is valid to all orders of QCD perturbation theory. Furthermore this result remains valid 
for all values of hadronic center-of-mass rapidity, because we have introduced a suitable integral 
transform over rapidity. The resummed rapidity-integrated Drell-Yan coefficient function to 
NLL is well known [U [2] and, using the notation of [lOj, it is given in a compact form (in the 
MS scheme) by 



Cr(iV,^,«V))=exp]-/ - l^-[A,a.i-)^A,aH-))^B,a.i':L 



Q'dk' ( , ,k\ , , _2,«:%\ „ ,Q 

nij? 



,(27) 
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where 



A, 



TT 



27r2 



67 



TT 
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2ti 



{21 



with Cp = 4/3, Ca = 3, Nj the number of flavors and with the Euler gamma = 0.5772 . . . . 
The use of only the first coefficient Ai allows us to resum all the LL contributions ag\og^~^^{N) 
while the use of all the three coefficients in eq. fl28|) enable us to add also the NLL terms 
«^log'(iV). 

A NLL expression of the rapidity distribution is obtained by taking the inverse Mellin and 
Fourier transform of a'^'^'^{N, Q^, M). This procedure requires the use of some specific prescription 
[T2| \T3\ [Ti] in order to overcome the problem of the Landau singularity in as{Q'^ /N"^)- Here, 
we adopt the "Minimal Prescription" proposed in [12], which is simply obtained choosing the 
integration contour of the inverse Mellin transform in such a way that all the poles of the 
integrand are to the left, except the Landau pole. Furthermore, in order to improve numerical 
convergence and to avoid the singularities of the parton densities of eq. (l26l) which are computed 
out of the real axis, we perform the A^-integral along a path F given by: 



F = F1 + F2 + F3 
M 

Fi(t) = Ca/p-^y +t(l + 0, te(-oo,0) 
M 

r2(s) = ^AfP + ^Sy, SG(-1,1) 

M 

TsW = CAfP + 2y-t(l-2), tG(0,+Oo) 



(29) 
(30) 

(31) 

(32) 



where Cmp is a positive number below the Landau pole of as{Q^ /N^). Performing the changes 
of variable M = — Inm and t = —Ins, the double inverse transform over the curve F becomes: 



(33) 



1 rl rjryi 

a''\x,Q^,Y) = - / cos(-Flnm)a""^(a;,g^-logm), 

TT JO m 

where (t'''^''(x, Q^, M) is given by 

a^^^(x,Q2,M)= (34) 

-Cmp -In s+i(M/2+l)^res(^^^^ ^ ^ _ •(^/2 + l),g^M)(l - i) 



1 

TT Jo 



1 ds 



+ 



sM 



-X 



-CMP~'i'SM/2 res 



a'-'\CMP + isM/2,Q\M) 



Eqs.f l33|IMIl are the expressions that we use to evaluate numerically the resummed dimension- 
less cross section in the variables x and Y . The explicit expression of eq. (l27I) is easily obtained 
performing the integrals and is given by 



N, — ,a^(/i^) 



exp{lnNgi{X)+g2{X)}, 



(35) 
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where 



[2A+(l-2A)log(l-2A)] 



(36) 



-'^-^ log(l - 2A) + ^[2A + log(l - 2A) + i log^(l - 2A)] 



-^[2A + log(l 
Po 



and 



A = /?oas(/i^)lnAr, 
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-Nf 
3 



167r2 



102 - -iV, 
3 ^ 



(37) 



(38) 



Here, we choose the factorization scale equal to the renormalization scale for simplicity. To 
study the dependence on the renormalization scale one has simply to express o;s(/x^) in terms of 
it. Furthermore, we need the analytic continuations to the whole complex plane of the Mellin- 
transformed parton densities that appear in eq. (l26|) . In order to overcome this problem, we have 
to evolve up a partonic fit taken at a certain scale solving the DGLAP evolution equations in 
Mellin space [I^. The LO and NLO expressions of the splitting functions are reported in [16] 
and their analytic continuations are given in [17] and |18j . 

Finally, we want to obtain a NLO determination of the cross section improved with NLL 
resummation. In order to do this, we must keep the resummed dimensionless part of the cross 
section eq. fl33l) . multiply it by the correct dimensional prefactors looking eqs. flimi"7|) . add the 
full NLO cross section and subtract the double-counted logarithmic enhanced contributions. 
This matching has to be done in the x and Y spaces, because we are not able to calculate the 
Mellin-Fourier moments of the full NLO cross section analytically. Thus, we have 



da 



da 



FO 



+ 



da' 



dQ'^dY dQ'^dY dQ^dY 



da' 



dQ^dY 



- a. 



d 



da' 



das KdQ'^dY 



(39) 



The first term is the full NLO cross section reported in [HI [13 [2D1 EI], which includes even the 
quark-gluon channel. The third and the fourth terms in eq. fl39p are obtained in the same way 
as the second one, but with the substitutions 



CT 



N, — ,a,(/i^) 



2^E - log 



(40) 
(41) 



respectively. The terms that appear in eq. fl41l) are exactly the 0{as) logarithmic enhanced 
contributions in the MS scheme. 

We note that the final expression eq. (l39|) is relevant even when the variable x is not large. In 
fact, the cross section can get the dominant contributions from the integral in eq. fllSI) for values 
of z eq.(l6]) that are near the threshold even when x is not close to one, because of the strong 
suppression of parton densities Fii^Xi,^"^) when Xi are large. 
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Figure 1: Y-dependence of (fa/{dQ'^dY) in units of pb/GcV^. The curves are, from top to bottom, the NLO 
result (red band), the LO+LL resummation (blue band) and the LO (black band). The bands are obtained 
varying the factorization scale between /i-^ = 2Q^ and = 1/2Q^. 



To show the importance of this resummation, we have calculated the Drell-Yan rapidity 
distribution for proton-proton collisions at the Fermilab fixed-target experiment E866/NuSea 
[22]. The center-of-mass energy has been fixed at a/S* = 38.76 GeV and the invariant mass 
of the virtual photon 7* has been chosen to be = 64 GeV^ in analogy with [TT] . Clearly 
the contribution of the virtual can be neglected, because its mass is much bigger than Q^. 
In this case x = 0.04260 and the upper and lower bound of the hadronic rapidity Y eq.([9]) 
are given by ±1.57795. We have evolved up the MRST 2001 parton distributions (taken at 
jj? = 1 GeV^) as in [llj, where the NNLO calculation is performed. However, results obtained 
using more modern parton sets should not be very different. The LO parton set is given in [23j 
with a^'^^mz) = 0.130 and the NLO set is given in [23] with a^^^{mz) = 0.119. The evolution 
of parton densities at the scale fi^ has been performed in the variable flavor number scheme. 
The quarks has been considered massless and, at the scale of the transition of the flavor number 
{Nf — >■ Nf + 1), the new flavor is generated dynamically. The resummation formula eq. (l26ll 
together with eqs. (135113811 has been used with the number of flavors Nf = 4. 

In figure[I], we plot the rapidity-dependence of the cross section at LO, NLO and LO improved 
with LL resummation. The effect of LL resummation is small compared to the effect of the full 
NLO correction. We see that, at leading order, the impact of the resummation is negligible 
in comparison to the NLO fixed-order correction. This means that the NLL resummation is 
necessary. 
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Figure 2: Y-dependence of (P a / [dQ"^ dY) in units of pb/GeV^. The curves are, from top to bottom, the NLO 
result (red band), the NLO+NLL resummation (green band) and the LO (black band). The bands are obtained 
as in figure [1] 




Figure 3: Y dependence of d'^cr/ (dQ'^dY) in units of pb/GeV^. The curves are, from top to bottom, the NLO 
result (red band) and the NLO+NLL resummation (green band) together with the E866/NuSea data. The bands 
are obtained as in figure [T] 
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The LO, the NLO and its NLL improvement cross sections are shown in figure [2J The effect 
of the NLL resummation in the central rapidity region is almost as large as the NLO correction, 
but it reduces the cross section instead of enhancing it for not large values of rapidity. Going 
from the LO result to the NLO with NLL resummation, we note a reduction of the dependence 
on the factorization scale i.e. a reduction of the theoretical error. It is interesting to observe that 
logarithmically enhanced and constant terms account for more than 80% of the NLO contribution 
for all relevant rapidities. Therefore, they have the same sign. Nevertheless a suppression arises 
due to the shift in the complex plane of the dominant contribution of the resummed exponent. 
This suppression starts at order O(a^). 

In figure [31 we report the experimental data of [22] converted to the Y variable, together 
with our NLO and NLL resummed predictions. The agreement with data is good and a great 
improvement for not large rapidity is obtained with respect to the NLO calculation. We note 
also that the NLL resummation gives better result than the NNLO calculation performed in [TT] . 
The NNLO prediction has a worse agreement with data than the NLO one for not large values 
of rapidity. This result suggests that, for the case of rapidity distributions, NLL resummation 
is more important than high-fixed-order calculation and that it can be so even at higher center- 
of-mass energies. 

To summarize, we have proved a resumation formula for the Drell-Yan rapidity distributions 
to all logarithmic accuracy and valid for all values of rapidity. Isolating a universal dimension- 
less coefficient function, which is exactly that ones of the Drell-Yan rapidity-integrated, we have 
shown a general procedure to obtain resummed results to NLL for the rapidity distributions of 
a virtual photon 7* or of a real vector boson W"^, . Furthermore, we have outlined a general 
method to calculate numerical predictions and analyzed the impact of resummation for the fixed- 
target experiment E866/NuSea. This shows that NLL resummation has an important effects on 
predictions of differential rapidity cross sections giving an agreement with data that is better 
than NNLO full calculations. 
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